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ABSTRACT 

We use observations of six "blazars" with the Very Long Baseline Array 
(VLB A), at eight frequencies (4.6, 5.1, 7.9, 8.9, 12.9, 15.4, 22.2, 43.1 GHz), 
to investigate the frequency-dependent position of their VLBI cores ("core- 
shift") and their overall jet spectral distribution. By cross-correlating the op- 
tically thin jet emission, we are able to accurately align the multi-frequency 
images of three of the jets (1418+546, 2007+777, 2200+420), whose core-shifts 
and spectra we find consistent with the equipartition regime of the Blandford 
& Konigl conical jet model, where the position of the radio core from the base 
of the jet follows r core (v) oc v~ x . For the jet of 0954+658, we align the higher 
frequency images using our lower frequency measurements assuming equipar- 
tition in the radio core from 4.6-43 GHz. The jet emission of the other two 
sources in our sample (1156+295, 1749+096) is too sparse for our alignment 
technique to work. Using our measured core-shifts, we calculate equipartition 
magnetic field strengths of the order of 10's to 100's of mG in the radio cores 
of these four AGN from 4.6-43 GHz. Extrapolating our results back to the 
accretion disk and black hole jet-launching distances, we find magnetic field 
strengths consistent with those expected from theoretical models of magneti- 
cally powered jets. 

Key words: radio continuum: galaxies - galaxies: jets - galaxies: magnetic 
fields 
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1 INTRODUCTION 



Highly collimated jets of relativistic plasma are considered to be generated and propelled 
outwards from the central regions of Active Galactic Nuc lei (AGN) by magnetic forces 



Meier 



2009, and references therein). High 



surrounding the supermassive black hole (see 
resolution VLBI observations can detect the synchrotron emission from these jets in the 
radio regime at distances of the order of 10 3 — 10 7 r g , where r g = GM BH /c 2 , correspond- 
i ng to the sub-parsec to parsec scales of the iet ([Lobanovi 120061 ). Current theoretical models 



(IVlahakis fc Konigl 



2004; 



Komissarov et al. 



20071 . hereafter K07) invoke rotation of the black 



hole magnetosphere and/or the inner accretion disk to create a stiff helical magnetic field 
that expels, accelerates and helps collimate the jet flow; an external medium is still re- 
quired to confine t he jet at its boundary, with a p ossible solution being a slower disk - driven 
win d outflow (e.g.. iBogovalov fc Tsinganosl 120051 ) . The models of 



and 



Blandford fc Reed (119741 ) 



Daly fc Marscherl (119881 ) show how jets can be accelerated hydrodynamically, however, 



Konigll (120071 ) proposes that extended acceleration regions are distinguishing characteris- 



tics of magnetohydro dynamic (MHD) dominate d flows over purely hydrodynamical ones, 



and uses the model of 



Vlahakis fc Konigl (|2004l ) to explain the parsec-scale acceleration of 



component C7 of 3C345 (jUnwin et al.lll997l ). 



In their conical jet model, iBlandford fc Konigll (119791 ) [hereafter BK79] proposed that 



"the unresolved core be identified with the innermost, optically thick region of the approach- 
ing jet" with the base of the jet corresponding to the vertex of the cone. The radio cores 
of AGN generally have flat or inverted (a > 0) spectra, where S cx v a . The optical depth 
due to synchrotron self-ab sorption (SSA) is described by r(r) oc r { 1 ^- a ) m + n - l h , 2 ^- a (e.g., 



Rybicki &: Lightmanlll979l ). where m and n describe the power-law fall-off in the magnetic 



field strength (B = Bi pc r~ m ) and particle number density (N = Ni pc r~ n ), respectively, 
with distance from the central engine. Since VLBI resolution is usually not sufficient to 
completely resolve the true optically thick radio core, the core region contains emission from 
around the r = 1 surface at that frequency, as well as some contribution from the optically 
thin inner jet. 

Furthermore, due to the frequ ency depende nce of the r = 1 surface, the position of 



the radio core follows r oc v x l kr ( IKonigl 



19811 ). where r is the distance from the central 



engine and k r = ((3 — 2a)m + 2n — 2)/(5 — 2a). For equipartition between the jet particle 
and magnetic field energy densities, the quantity k r = 1, with the choice of m = 1 and 
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n 



2 being physically reasonable (e.g.. lKonigl 



1981 



Hutter fc Mufson 



1986 



Lobanov 



1998a . 



hereafter L98), making k r independent of the spectral index. Therefore, in equipartition, the 
"core-shift" (Ar = r Ul — r U2 ) between two frequencies {y-i > Ui) is directly proportional to 
(z/ 2 - v x )j{y 2 v x ). 

Since absolute positional information is lost during calibration of VLBI images, they 
are centered at the peak intensity of the map, which does not necessarily correspond to 
the furthest upstream component referred to as the VLBI core. For accurate combination 
of multi-frequency VLBI maps, we must properly align our images by correcting for the 
frequency-dependent posit ion of the self-absorbed radio core. Phase-reference observations 
( IMarcaide fc Shapiro! 1 1984 ) can be used to obtain the absolute position of the radio core by 
comparing the position of the target source with respect to a reference source, which should 
ideally be a point source, but this is a very calibration-intensive and weather dependent 
approach and gets much more difficult with observations at more than one frequency. The 
frequency dependent core-shift has generally been obtained by aligning model-fitted opti 
cally thin jet comp onents at different frequencies ("self-referencing") (e.g., 



Lobanov 



1998a 



Kovalev et al. 



20081 ). but this method is unsuitable for smooth or weak jets that do not have 



a distinct component across a wide range of frequencies on which to align the images. An 
alternative method of aligning VLBI images is to cross-correlate the optically thin emission 
from the jet at multiple pair s of frequencies from regions where we do not expect large vari- 



ations in the spectral index. 



Walker et al 



(120001 ) used such a cross-correlation technique to 
successfully align his multi-frequency images of 3C84. We also use a cross-correlation tech- 
nique for aligning our images, which involves comparing the optically thin jet emission at 
two frequencies with the highest co r relatio n-coefficient occurring when the two images are 



best aligned; see 



Croke fc Gabuzda 



(120081 ) for a detailed description of the method. From 



this method, we have been able to find reliable measurements for t he frequency dependent 



core- shift of several sources across a wide range of frequencies (see 



P'Sullivan fc Gabuzda 



2008 



, for comparisons with both phase refere ncing and "self- r eferencing" methods). 
From long-term monitoring of blazar jets, iMarscherl (120091 ) lists evidence supporting the 
i dea that in some cases the cores correspond to a stationary feature such as a conical shock 
( 1Cawthornell2006l ) . He states that the core would correspond to either the r = 1 surface or the 
first stationary hotspot downstream of that surface, whichever is more intense at a particular 
epoch. Our core-shift measurements after the alignment of the multi-frequency images of 
2200+420, 2007+777, 1418+546 and 0954+658 are consistent with the assumption of the 
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Table 1. Model Parameters for all sources 



Source 


Z 






<$> 


r 


e 


8 






(Mpc) 




(°) 




(°) 




0954+658 


0.368 


1958 


-0.6 


2.2 


6.3* 


9.3 


6.2* 


1156+295 


0.729 


4481 






25.1 1 " 


2.0* 


28.5* 


1418+546 


0.152 


716 


-0.5 


1.9 


4.5* 


11.2 


5.1* 


1749+096 


0.320 


1662 






7.5t 


3.8* 


12.0* 


2007+777 


0.342 


1797 


-0.5 


1.0 


8.0 


7.2 


7.9* 


2200+420 


0.069 


307 


-0.5 


1.9* 


7.0* 


7.7* 


7.2* 


Cosmology: 


H = 


71 km s — 1 


Mpc" 




= 0.73 


and Qm 


= 0.27. 



luminosity distance, Col. (4): jet spectral index, Col. (5): jet half-op e ning an gle, Col. (6): bulk Lore n tz fac to r, Col. (7): jet 
viewin g angle, Col. (8 ) : Dop pler factor. References: iGabuzda et al.1 lll994T) . jPushkarev fc Gabuz di Jl999h . iHovatta et al.l 
l|2009h . $ Jorstad et at] (|2005h 



core corresponding to the r = 1 surface and not some stationary feature. It's also worthwhile 
inspecting the core polarization because this can sometimes give information about scales 
smaller than the standard re solution at a particular f reque ncy. From the polarization images 



of these sources presented in lO'Sullivan fc Gabuzda (120091 ) , we do not see any clear evidence 



of distinct polarized features upstream of the total intensity core which would support the 
model of a stationary core feature through which ejected plasma blobs propagated. 

In Section 2 we describe our observations and how we use our core-shift measurements 
to calculate magnetic field strengths. Our results for each source are presented in Section 3 
and we discuss some of their interesting features in Section 4. Our conclusions are listed in 
Section 5. 



2 OBSERVATIONS AND METHOD 

Six "blazars" were observed on 2006 July 2 with the American Very Long Baseline Ar- 
ray (VLBA) at ei ght frequencies (4.6, 5.1, 7. 9, 8.9, 12.9, 15.4, 22.2, 43.1 GHz) quasi- 



simultaneously (see lO'Sullivan &: Gabuzdall2009l . for details). Here we present spectral index 



images for these six blazars along with estimates of the magnetic field strengths in the radio 
core from 4.6-43 GHz for four sources in which we were able to measure the "core-shift" . We 
follow the method described in L98 for calculating several important jet parameters from 
the frequency-dependent shift of the position of the VLBI core due to SSA. Tabled] lists the 
observed sources along with the model parameters used in our calculations. 

Once our images are properly aligned, we generate spectral index (a) maps between 
pairs of frequencies in aips using the task comb, with a = (log S U2 — log S Ul )/ (log z/ 2 — log V\) 
where S u is the measured flux at frequency v. This approximates the spectrum as a straight 
line which should be a good approximation as long as there is not a spectral turnover 



Magnetic Field Strength in AGN Jets 5 

betw een the two freq uencies. Spectral turnovers of components are often observed in VLBI 
jets (jLobanovlll998bl ) but useful information can still be extracted by generating maps from 
relatively closely spaced frequency pairs. In order to analyse the spectrum across our wide 
frequency range, we constructed spectral index maps from 4.6-7.9, 8.9-12.9, 15.4-22 and 
22-43 GHz for each source. The images used to create each map have the same image size, 
pixel size and are convolved with the same beam, corresponding to the lower frequency of 
the pair. We also take a slice along the jet of each spectral index image for inspection. 

Since the shift measurements from the cross-correlation program are found with respect 
to the Stokes / peak intensities, we use them in conjunction with the model-fitted core 
positions from difmap ( IShepherdlll9971 ) to obtain the distance between the cores at different 
frequencies ("core-shift"). The task uvrange in difmap was used to provide nv-ranges for 
the higher frequencies matching that of the lower frequencies when model-fitting. Due to 
the fall-off in the optically thin part of the synchrotron spectrum with increasing frequency 
and the large differences in resolution, aligning images across wide spacings in frequency can 
prove very difficult. However, with our excellent frequency coverage from 4.6-43 GHz, we 
are able to obtain many measurements from different frequency combinations which can be 
compared with each other to check for consistency. 

To calculate the ma gnetic field stren gth with distance from the base of the jet, we combine 
the analysis of L98 and iHirotanil (120051 ). The distance of the radio core from the base of the 
jet at the observed frequency v in GHz, is given by 

9^) v ~l/kr (1) 

smO J 

in parsecs, where 6 is the jet viewing angle and Q ru is the core-position offset measure, 
defined in L98 as 

O - 4 85 x IP' 9 Ar ™° DL ( vl /kr »2 /kr \ (2) 
Or, -4.85x10 {1 + z)2 y u x /kr _ u i/k r ) ( 2 ) 

in units of pc-GHz, where is the luminosity distance in parsecs and Ar mas is the measured 
core-shift in milliarcseconds. 

We use the synchrotron self-absorption derivation for a conical jet by 
and the equipartition condition, in cgs units, 

B 2 

N^ min m e c 2 = K— (3) 

<37T 

K _ 2a + l (Tmax/Tmin) 2 " ~ 1 ( ) 

2a ( 7 max/7min) 2a+1 - 1' 1 J 

where 7 m i n and 7 max are the low and high cut-offs of the electron energy distribution. Eval- 



Hirotanil (120051 ) 
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uating eqn. 43 from iHirotanil (120051 ) for the simplified case k r = 1 and a 
equipartition magnetic field in Gauss at 1 pc 

'^(1 + *) 2XV4 



B x ~ 0.025 



sin 2 fl 



-0.5, we get the 



(5) 



where 5 is the Doppler factor, 9 is the viewing angle and (f) is the half-opening angle of the 
jet. The magnetic field strength in the core at the observed frequency v can then be found 
from 



B r 



^1' core 



(6) 



For spectral indices other than —0.5, the magnetic field strength depends on the exact 
value of 7 m i n according to B\ oc r )^^ a ^ <y7 ^ • The dependence on K is much weaker since 
7max 3> 7min- Equation [5] differs from the equipartition magnetic field at 1 pc derived in 
L98, eqn. 10, where our calculated magnetic field strengths are a factor of ~ 3 times smalle r 
(e.g., compare our results for 2200+420 here with those in 
due to the more accurate equipartition condition used in 



O 'Sullivan k Gabuzdal (|2008j)) 



Hirotani 



(120051 ) . Note that in 



all cases, where errors are presented, we have only considered the errors in the core-shift 
measurements. 



3 MAGNETIC FIELD STRENGTH AND SPECTRAL INDEX RESULTS 
3.1 2200+420 (BL Lac) 

Fig. [1] shows the spectral index distribution for 2200+420 from 4.6-43 GHz. The core region 
displays a strongly inversted spectrum at low frequencies (Fig. [T^l) before becoming less 
inverted at higher frequencies (Fig. [lb, c, d), consistent with SSA of emission from the 
unresolved part of the jet. Further downstream, the jet emission becomes optically thin, 
with the majority of slices consistent with a spectral index of ~ —0.5. The direction of 
the slice along the core region in the 22-43 GHz map (Fig. [TH) follows the 43-GHz inner jet 
direction. There is an interesting feature in the 8.9-12.9 GHz map (Fig. [lb), at ~5 mas from 
the core, where the spectrum flattens considerably. This feature is also present in the 4.6-7.9 
GHz and 15.4-22GHz maps, but is less pronounced, possibly due to the low resolution at 4.6 
GHz and the lack of sensitivity to the extended jet emission at the higher frequencies. The 
presence of this feature might be due to some particle re-acceleration in the jet or interaction 
between the jet and the surrounding medium, since this is close to where the jet changes to 
a south-easterly direction in the plane of the sky. 
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The measured core-shifts for 2200+420 are listed in Table [2] along with the angle repre- 
senting the direction of the shift measured from north through east. The average direction 
of the shift is 14 ± 10°, where the error corresponds to the standard deviation, and this 
core-shift direction is consistent with the inner jet direction at 43 GHz of ~ —163°. Table [3] 
lists the size of the shift as a fraction of the beamsize. Even though the shifts are only of 
the order of 5-20% of the corresponding beams, they still have quite a dramatic effect on 
the spectral index in the core if they are not taken into account. 

Figure [2] shows each shift measurement, Ar, plotted against the expected frequency 
dependence for k r = 1. The y-intercept from a straight line fit to the data is —0.01 ± 0.01 
mas. From the slope of the line we determine a characteristic core position offset measure 
(fl ru ) of 3.4 ±0.3 pc ; GHz for 2200+420 . Our results are also consistent with a previous core- 



shift estimate from 



Mutel et al 



Jl990j) who found a shift of 0.3 mas between 5.0 and 10.6 
GHz, corresponding to Q ru = 3.7 pc • GHz (no error was presented with their result). Since 
the position of the core is expected to change with frequency according to r corc (z/) oc z/~ 1//fcr , 
observations at three or more frequencies enable determination of the parameter k r . Our 
dense frequency sampling allows us to provide strong constraints on the value of k r across 
our entire frequency range. Using 43 GHz as the reference frequency and averaging the 
various possible combinations of shift measurements (Table [3]), we find our values to be 
consistent with a constant k r of 0.99 ± 0.07 from 4.6-43GHz (Fig. [3]). This indicates that 
the physical conditions in the compact jet region in 2200+420 do not change significantly 
between 4.6 and 43 GHz. A value of k r = 1 indicates that the jet particle and magnetic 
field energy densities are in equipartition, and our derived value of k r is consistent with this 
assumption. 

Using our core-shift measurements, we can find estimates for the magnetic field strength, 
distance of the radio core to the base of the jet and the total jet luminosity. To do this 
we follow the method outlined in L98. Estimates of the jet opening angle, viewing angle, 



bulk Lorentz 
2200+420 by 



actor and Doppler factor are taken from multi-epoch VLBI observations of 
(200 5,1), see Table [D These estimates are also consistent with 



Jorstad et al. 



more recent observations by 



Hovatta et al. 



(120091 ) , who analyse the variability of radio flares 
from single dish measurements to obtain their jet parameters. We use a jet spectral index 
(ctj) of —0.5 for our calculations, consistent with our spectral index maps for this source 
(Fig. [T|). We use the shift measurements listed in Table [2] to calculate the average core- 
position offset measure (fl r „) from eqn. ([2]) at each frequency. Then for each Q ru , we calculate 
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Table 2. All measured core-shifts for 2200+420, plotted in Fig.|2] 



^1 


^2 


Ar 


Angle 


(GHz) 


(GHz) 


(mas) 


(°) 


4 6 


5 1 


n nfi -I- n n7 

u.uu in U.U 1 


— 1.0 zb 0.1 


4 6 


7.9 


n i q -t- n 07 

U. -L y HI KJ.U 1 


J-O. ^ zc u. o 




8 9 


n on 4- n (17 
u. ^u in U.U/ 


1 r: c _L n o 
J-O. u zc u. o 


4 5 


12.9 


n or: _|_ n fl7 
u. oo zn U.U/ 


1Q4 + 04 


4.6 


15.4 


40 + 07 

VJ.'-tVJ _1_ VJ.KJi 


i k q J- n q 


5 1 


7.9 


n 1 4 -i- n 07 

U. it ir_ u.u/ 


.^z . u zc u.u 


5 1 


S 9 


u. j. y ir_ u.u/ 


16.4 zb 0.3 


5 1 


12.9 


n Qfi _|_ n n7 
u.ou ir_ u.u/ 


z,o.Q in u.u 


5 1 


15.4 


U.OJ ir_ u.u/ 


19.1 i 0.4 


7.9 


S 9 


U.U'i IT_ U.U'i 


5.6 i 0.1 


7.9 


12.9 




±o. o zc u. o 


7.9 


15.4 


0.18 ±0.04 


14.0 ±0.1 


7.9 


22.2 






7.9 


43.1 


0.28 ±0.04 


6.4 ±0.1 


8.9 


12.9 


0.07 ±0.04 


11.0 ±0.1 


8.9 


15.4 


0.11 ±0.04 


20.9 ±0.3 


8.9 


22.2 


0.17 ±0.04 


13.1 ±0.1 


8.9 


43.1 


0.18 ±0.04 


-0.6 ±0.1 


12.9 


15.4 


0.03 ± 0.03 


-9.1 ±0.1 


12.9 


22.2 


0.07 ±0.03 


12.7 ±0.1 


12.9 


43.1 


0.13 ±0.03 


2.3 ± 0.1 


15.4 


22.2 


0.05 ± 0.03 


23.5 ±0.3 


15.4 


43.1 


0.09 ± 0.03 


19.9 ±0.2 


22.2 


43.1 


0.03 ± 0.02 


33.2 ±0.4 



Col. (1) & (2): observed frequency in GHz, Col. (3): magnitude of the core-shift in milliarcseconds, Col. (4): angle of core-shift 
in degrees. Average angle: 14 ± 10° 



Table 3. Averaged core-shifts using 43 GHz as the reference frequency for 2200+420, plotted in Fig. [3] 



V 


Ar 


A?" projected 


Fraction of beam 


(GHz) 


(mas) 


(pc) 


(%) 


4.6 


0.46 ± 0.07 


0.60 ±0.09 


17 


5.1 


0.43 ± 0.07 


0.56 ±0.09 


17 


7.9 


0.27 ±0.04 


0.35 ±0.05 


16 


8.9 


0.21 ±0.04 


0.27 ±0.05 


14 


12.9 


0.12 ±0.03 


0.16 ±0.04 


11 


15.4 


0.10 ±0.03 


0.13 ±0.04 


10 


22.2 


0.04 ±0.02 


0.05 ±0.03 


5 



Col. (1): observed frequency in GHz, Col. (2): magnitude of the shift from the 43-GHz core position, Col. (3): projected shift 
in parsecs (1 mas = 1.3 pc), Col. (4): core-shift as a fraction of the beam-size at the corresponding frequency, in percent. 



the de-projected distance of the radio core to the base of the jet, r core {y) from eqn. ([T]), the 
equipartition magnetic field strength at 1 pc, B\ from eqn. ([5]), and finally, using eqn. ([6]) 
the magnetic field strength in the radio core. Our results are listed in Table HI 

It's important to note the effect that small changes in the parameter k r have on our 
calculated values. We get 0.92 < k r < 1.06 from the range of possible values within the 
error in the fit shown in Fig. [3j There is a strong dependence on the value of r core [y) with 
3.5 < r core (43 GHz) < 0.14 pc, while the magnetic field strength at 1 pc remains relatively 
insensitive to small variations in k r with 0.12 < B\ < 0.14 pc (using Q ru at 43 GHz). 
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Therefore, while there is a rather large uncertainty in r core and hence B core , we consider the 
calculated value of B\ to be reliable for the jet parameters listed in Table [U 

Adopting k r = 1, we plot the magnetic field strength in the core as a function of the 
distance from the base of the jet in Fig. HI Fitting the line B(r) = B\r~ l to the data, we 
find B\ = 0.139 G with a formal error of 0.006 G. By extrapolation, our results predict 



field strengths of the orde r of 1 G in the 



results for other jets from 



Savolainen et al. 



mm core of 2200+420. This is consistent with 



(120081 ). where he finds magnetic field strengths 



of order 1 G in the core of 3C 273 by measuring the synchrotron self- absorption turnover 



frequency (e.g. 



Marscher 



19871 ) of individual model-fitted components from 5-86 GHz (60- 



3.5 mm). Importantly, this method does not require the assumption of equipartition to 
separate the magnetic field strength from the particle number density, but the calculations 
are very sensitive to the exact values of the turnover frequency and the estimated size of the 
model-fitted core component at the turnover frequency. 

In their simulations of magnetically driven relativistic jets, K07 consider applications of 
their model to AGN jets and provide a means of estimating the distance from the central 
engine at which equipartition between the Poynting and matter energy fluxes is established, 
-Req,BH ~ 2x 10 16 (Mbh/1O 8 M ) (0. 1/0) cm, for jet s launched from the black hole or from the 
nuclear accretion disk R eq (AD) ~ 1O 17 (Mbh/1O 8 M )(O.1/0) cm. In their model, the jet flow 
is magnetically accelerated up to a distance of ~ 10R eq , at which point it reaches the terminal 
Lorentz factor; hence, observations that probe this region may detect deviations from the jet 
model of BK79, which assumes a constant jet speed, and possibly infer whether the observed 
jet has been launched from the black hole or the accretion disk. Using the jet opening angle 



of 1.9° and a black hole mass of 1.7 x 1O 8 M (I Woo fc Urrv 



2002J), the equipartition distance 



for a black-hole-launched jet in 2200+420 is ~ 0.03 pc while R eq ~ 0.2 pc for a jet launched 
from the accretion disk, which according to our results, would be accessible with VLBI 
observations at frequencies greater than 120 GHz. This is out of the range of cur rent VLBI 



Doeleman et al. 



2009J). However, 



but could be probed with future sub-mm VLBI arrays (e.g., 
since the terminal Lorentz factor is not reached until ~ 10R eq it's possible that accelerating 
jet components could be seen in the inner jet at frequenci es as low as ~ 15 GHz. Multi-epoch 



VLBA observations at 43 GHz fromlJorstad et al. 



(120051 ) have indeed found accelerating jet 



components for this source which they attribute to a small change in the Lorentz factor 
along with a substantial change in the direction of the jet flow. The large uncertainty in 
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Table 4. Derived physical quantities for 2200+420 (k r = 1) 



V Bcore ) TcQT&{y) -^lpc £lrv 

(GHz) (G) (pc) (G) (pc-GHz) 



4.6 


0.03 ±0.01 


5.9 ±5.0 


0.16 ±0.20 




5.1 


0.03 ±0.01 


5.4 ±4.5 


0.16 ± 0.20 


3.7 ±3.1 


7.9 


0.05 ±0.01 


2.5 ±0.7 


0.13 ± 0.05 


2.8 ±0.7 


8.9 


0.06 ±0.01 


2.5 ±0.6 


0.14 ±0.05 


3.0 ±0.8 


12.9 


0.08 ±0.01 


1.9 ±0.4 


0.15 ± 0.05 


3.3±0.7 


15.4 


0.09 ±0.01 


1.6 ±0.3 


0.15 ± 0.04 


3.4 ±0.6 


22.2 


0.14 ±0.02 


1.1 ±0.3 


0.14 ±0.06 


3.1 ±0.9 


43.1 


0.27 ±0.02 


0.5 ±0.1 


0.13 ± 0.03 


2.8 ±0.4 



Col. (1): observed frequency in GHz, Col. (2): radio core magnetic field strength in Gauss, Col. (3): de-projected distance of 
the radio core to the base of the jet in parsecs, Col. (4): magnetic field strength at 1 pc, Col. (5): average core-position offset 
measure, in units of pc ■ GHz, for shifts with V2 hi Table [2] equal to v in Col. (1) above. Note: Physical parameters at 4.6 GHz 
calculated with 5.1 GHz core-position offset measure. 

r core should be kept in mind though, as r core (43 GHz) may extend as far as 3.5 pc from the 
base of the jet. 

Using the divergence theorem, we can also estimate the initial magnetic field strength 
at the jet launching distance by assuming the energy extracted per time is equal to the 
Poynting flux times the area (N. Vlahakis, private communication), which gives B ~ 1.5 x 
lO- 5 /, 1 ^- 1 . For example, with L ~ 10 46 erg s" 1 and M BH ~ 10 8 M , we get B ,ad ~ 10 4 G 
at an inner accretion disk distance of 10r 9 and Bq(BH) ~ 10 5 G at a black hole jet-launching 
distance of r g , as quoted in K07. 

Using eqn. (23) from BK79, we can find an estimate for the total power carried by the 
relativistic electrons and magnetic field in the jet, Ltotai ~ 2.5 x 10 A8 T 2 /3(f) 2 B 2 erg s _1 , where 
(3 is the speed of the jet in units of c. From our best-fit value of B\ ~ 0.14 G, we get 
^totai ~ (2.6 ± 0.2) x 10 45 erg s _1 . Hence, for 2200+420, we expect B (AD) ~ 3 x 10 3 G 
and Bq(BH) ~3x 10 4 G. In Fig. we extend our plot of B vs. r to include the calculated 
values of the initial magnetic field strength at the corresponding jet launching distances. 
Extrapolating 0.14r _1 all the way down to lr g , we find B ~ 1.7 x 10 4 G, in good agreement 
with the estimated magnetic strengths at the nominal jet launching distances. 

3.2 2007+777 

In Fig. [6] we show the spectral index distribution of 2007+777 from 4.6-43 GHz. For this 
source we find a flat to inverted spectrum in the core region, falling off to a steep optically 
thin spectrum further downstream in the jet. The slices for the 4.6-7.9 GHz and 8.9-12.9 
GHz maps (Figj6^, b) show optically thick regions at the edge of the jet ~8 mas from the 
core; however, these are possibly spurious features due to edge effects when combining the 
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two maps. We only show the spectral index map from 15.4-43 GHz due to loss of flux in the 
22-GHz map. 

All measured core-shifts for 2007+777 are listed in Table 5. The average direction of 
the measured shift is 89 ± 13°, which is consistent with the parsec scale jet direction of 
~ -90°. The y-intercept of the best-fit line to the data (Fig. ED is -0.001 ±0.004, indicating 
that the compact jet region from 4.6-43 GHz is close to equipartition. Due to a failure to 
find shift measurements that produced reasonable spectral index maps combining the higher 
frequencies with the 7.9 and 8.9 GHz images, we do not have as many measurements as in 
the previous case of 2200+420. In order to make the spectral index map between 8.9 and 12.9 
GHz, we use in Fig. [6]the characteristic position-offset measure of 5.2 + 0.3 pc-GHz from the 
best-fit line to the measured data in Fig. [7] to estimate the required shift. This produces a 
spectal index map consistent with those generated at lower and higher frequencies (Fig. [6]). 
In order to find an estimate for the parameter k r , we use over-lapping shifts to plot 7.9-43 
GHz and 8.9-43 GHz values (e.g., (4.6-15.4 GHz)-(4.6-7.9 GHz) + (15.4-43 GHz)). Hence, 
we find k r = 0.97 ± 0.06 (Fig. [8]) for this source (see Table 6 for listed values). Without the 
7.9 -43 GHz and 8 9-43 G Hz points, we get k r = 0.99 ± 0.05. 

Gabuzda et al.l (119941 ) find /3 app = 2.9 for 2007+777, where /3 app is the apparent velocity 



of the jet in units of the speed of light, based on three observations between 1980 and 
1989. Howe ver, from much mo r e den sely sampled multi-epoch monitoring from 1994-2002 



at 15 GHz, 



2007+777. 



Kellermann et al 



Hovatta et al 



(120041 ) do not detect any outward motion from the jet of 
(120091 ) derive a Doppler factor (5) of 7.9 using variability arguments 
from single dish monitoring indicating that this jet does have a relativistic flow. To estimate 
the magnetic field strength in this source, we require estimates for the bulk Lorentz factor 
(T) and the viewing angle (8). Due to the uncertainty in the jet speed measurements, we 



assume (3 ; 



app 



5 and from this we find the minimum Lorentz factor of the flow from r m j r 



(l+/3f p ) 1 / 2 . To estimate the viewing angle we use 9 = sin 1 (l/r min ); while this in unlikely to 
be tr ue for any one source, it is a good approximation in a statistical sense (see 



Cohen et al. 



20071 ). We find a jet opening angle of 1.0° by estimating the projected half opening angle 



using the ratio between the apparent transverse size and longitudinal distance of the jet 
components within 1 mas from the cor e (i.e., we assume a conical jet with a constant 
opening angle within this region, see also lJorstad et al.l (120051 )). 

Due to the absence of direct jet speed measurements, the following magnetic field values 
can only be considered as rough estimates. Using the same procedure as described in the 
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Table 5. All measured core-shifts for 2007+777, plotted in Fig.0 



vi V2 Angle 

(GHz) (GHz) (mas) (°) 



4.6 


5.1 


0.02 ±0.02 


91 ± 1 


1.6 


7.9 


0.10 ± 0.02 


89 ± 1 


4.6 


8.9 


0.12 ±0.02 


89 ± 1 


4.6 


12.9 


0.15 ±0.02 


79 ±1 


1.6 


15.4 


0.17 ±0.02 


87 ±1 


5.1 


7.9 


0.08 ±0.02 


97 ± 1 


5.1 


8.9 


0.10 ±0.02 


88 ± 1 


5.1 


12.9 


0.12 ±0.02 


82 ± 1 


5.1 


15.4 


0.13 ±0.02 


92 ± 1 


7.9 


8.9 


0.02 ±0.01 


91.2 ± 0.8 


12.9 


15.4 


0.016 ±0.007 


117.1 ± 0.3 


12.9 


22.2 


0.030 ± 0.007 


104.1 ± 0.4 


12.9 


43.1 


0.055 ±0.007 


84.4 ±0.4 


15.4 


22.2 


0.018 ±0.007 


54.8 ± 0.3 


15.4 


43.1 


0.045 ± 0.007 


86.7 ±0.4 


22.2 


43.1 


0.022 ± 0.006 


83.5 ± 0.3 



Col. (1) & (2): observed frequency in GHz, Col. (3): magnitude of the core-shift in milliarcseconds, Col. (4): angle of core-shift 
in degrees. Average angle: 89 ± 13° 



case of 2200+420, the calculated B core , r core and Bi are presented in Table [7] using the 
equipartition value of k r = 1 and a jet spectral index of —0.5 consistent with the spectral 
index maps (Fig. [6]). Considering a range of values for k r , from the fit to the data shown in 
Fig. El of 0.91 < k r < 1.03, we get 8.7 < r core (43 GHz) < 0.5 pc and 0.23 < B x < 0.27 pc 
(using Q ru at 43 GHz). 

We plot B core versus r core , for k r = 1, in Fig. [9] and find B\ = 0.257 ±0.007 G. Extending 
this line to the launching region of the jet (Fig. fTOl) . we find that it reaches ~ 1.3 x 10 3 G 
at an accreti on disk distance of 10r 9 and ~ 1.3 x 10 4 G at r g , using a black hole mass of 
4.0 x 10 8 M Q Jwu et al.lhoOQl ). We calculate the total jet power using B 1 ~ 0.26 G and with 
the jet parameters listed in Tabled] to find Ltotai ~ (3.2 ± 0.2) x 10 45 erg s _1 . Using L total , 
we estimate the magnetic field strength of B ~ 1.4 x 10 3 G at the jet launching distance of 
10r 9 , which is in excellent agreement with our extrapolated magnetic field strength results 

(Fig. no}. 

The estimated equipartition radius for an accretion-disk-launched jet of ~ 0.7 pc is just 
outside the range of 43-GHz VLBI, which probes distances of ~ 0.9 pc from the central 
engine (Fig. [TO]) . However, this source should be a good candidate for the possible detection 
of accelerating jet components from multi-epoch monitoring at 43 GHz, even though the 
43-GHz core may extend as far as ~ 9 pc from the base of the jet, since the jet flow is not 
expected to reach its terminal Lorentz factor until ~ 1 pc for a black-hole-launched jet and 
~ 7 pc for an accretion-disk-launched jet. 
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Table 6. Averaged core-shifts using 43 GHz as the reference frequency for 2007+777, plotted in Fig. [8] 



v Ar Ar pro j ectel i Fraction of beam 

(GHz) (mas) (pc) (%) 



4.6 


0.21 ±0.02 


1.0 ±0.1 


9 


5.1 


0.17 ±0.02 


0.8 ±0.1 


8 


7.9 


0.10 ±0.02 


0.5 ±0.1 


7 


8.9 


0.08 ± 0.02 


0.4 ±0.1 


6 


12.9 


0.056 ± 0.007 


0.27 ±0.03 


6 


15.4 


0.042 ± 0.007 


0.20 ±0.03 


6 


22.2 


0.025 ± 0.006 


0.12 ±0.03 


5 



Col. (1): observed frequency in GHz, Col. (2): magnitude of the shift from the 43-GHz core position, Col. (3): projected shift 
in parsecs (1 mas = 4.8 pc), Col. (4): core-shift as a fraction of the beam-size at the corresponding frequency, in percent. 



Table 7. Derived physical quantities for 2007±777 (k r = 1) 



Score (t^) Tcore (^) 



(GHz) 



(G) 



(PC) 



(G) (pc ■ GHz) 



1.6 


0.030 ±0.010 


8.7 ±4.9 


0.26 ±0.22 




5.1 


0.034 ±0.010 


7.8 ±4.4 


0.26 ±0.22 


5.0 ±2.8 


7.9 


0.051 ±0.003 


5.5 ±0.7 


0.28 ±0.05 


5.5 ± 0.7 


8.9 


0.056 ±0.003 


5.4 ±0.5 


0.30 ±0.05 


6.0 ±0.6 


12.9 


0.085 ± 0.005 


3.1 ±0.3 


0.26 ±0.04 


5.0 ±0.5 


15.4 


0.099 ± 0.005 


2.8 ±0.3 


0.28 ±0.04 


5.5 ± 0.5 


22.2 


0.15 ±0.01 


1.6 ±0.3 


0.24 ±0.06 


4.4 ±0.7 


43.1 


0.28 ±0.03 


0.9 ±0.2 


0.26 ±0.07 


5.0 ±0.9 



Col. (1): observed frequency in GHz, Col. (2): radio core magnetic field strength in Gauss, Col. (3): de-projected distance of 
the radio core to the base of the jet in parsecs, Col. (4): magnetic field strength at 1 pc, Col. (5): average core-position offset 
measure, in units of 10 9 pcGHz, for shifts with U2 in Table[5]equal to v in Col. (1) above. Note: Physical parameters at 4.6 GHz 
calculated with 5.1 GHz core-position offset measure. 



3.3 1418+546 



Fig. [Til shows the spectral index distribution from 4.6-43 GHz for the source 1418+546. We 
find an inverted core spectrum between the lower frequency pairs (Fig. PTTk. b, c) with a flat 
spectrum core region between 22 and 43 GHz (Fig. fTTti). In the 4.6-7.9 GHz map (Fig. [TTa) 
the jet spectral index is ~-0.2, while at higher frequencies it becomes steeper with a typical 
value of -0.5 (Fig. [Tib, c). 

The best fit line to the measured shifts listed in Table 8 gives a y-intercept of 0.003±0.004 
assuming k r = 1 and the slope provides a characteristic Q ru = 3.9 ± 0.2 pc ■ GHz (Fig. 12). 
The average position angle of the shift is —52 ± 13°, consistent within the errors with the 
inner jet direction at 43 GHz of ~ 137°. Again averaging the measured shifts and using 
43 GHz as the reference frequency (Table 9), we find k r = 1.08 ±0.11 between 4.6 and 43 
GHz (Fig. [13]). The equipartition value k r = 1 and ctj = —0.5 are used to calculate the 
physical parameters listed in Table 10. We estimate a jet opening angle of 1.9° from the 
size of the jet c ompo nents within 1 mas fro m the core. We use a Dqppler factor of 5.1 from 



Hovatta et al. 



(120091 ) and /3 app = 4.3 from iPushkarev fc Gabuzdal (119991 ) to find estimates 
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Table 8. All measured core-shifts for 1418+546, plotted in Fig. [T2l 



^1 


^2 


Ar 


Angle 


(GHz) 


(GHz) 


(mas) 


(°) 


4 6 


5 1 


U.UO ZZ KJ.Vl^t 


_QQ 9 4- Q 
oo.i zz u. y 


4 6 


7.9 


14 -1- fl 09 

KJ . ±"± ZZ KJ.Vl^t 


— 4Q fi+(17 


4 6 


8 9 


n i ^ 4- n 09 


—44 4- fi 


4 6 


12.9 


Li.iU zz u.u^ 


—51.7 ± 0.7 


1.6 


15.4 


0.23 ±0.02 


-60.3 ±0.5 


5.1 


7.9 


0.11 ±0.02 


-41.5 ±0.8 


5.1 


8.9 


0.13 ±0.02 


-35.8 ±0.6 


5.1 


12.9 


0.17 ±0.02 


-46.8 ±0.4 


5.1 


15.4 


0.21 ±0.02 


-54.0 ±0.5 


7.9 


8.9 


0.02 ±0.02 


-24.3 ±0.7 


7.9 


12.9 


0.07 ±0.02 


-53.2 ±0.7 


7.9 


15.4 


0.11 ±0.02 


-57.4 ±0.2 


8.9 


12.9 


0.06 ±0.02 


-58.4 ±0.9 


8.9 


15.4 


0.09 ±0.02 


-68.3 ±0.7 


12.9 


15.4 


0.021 ±0.007 


-82.4 ±0.3 


15.4 


22.2 


0.030 ± 0.007 


-62.7 ±0.3 


22.2 


43.1 


0.036 ± 0.005 


-56.2 ±0.2 



Col. (1) & (2): observed frequency in GHz, Col. (3): magnitude of the core-shift in milliarcseconds, Col. (4): angle of core-shift 
in degrees. Average angle: —52 ± 13° 

of the bulk Lorentz factor of 4.5 using; T = {j3l + 5 2 + 1)/(2S) and the je t viewing angle of 



Hovatta et al. 



20091 ). 



11.2° using 9 = tan- 1 ((2/3 app )/((/? a 2 pp + 6 2 - 1)) (e.g., 

We find B\ = 0.178 ±0.006 G from our plot of B core versus r core (Fig. [14]). Extrapolating 
this line, we find B ~ 680 G at an accr etion disk distance of 10r s , using the black hole mass 
estimate of 10 8 ' 74 from 



Wuetal 



(120091 ). Using the calculated total luminosity of (1.7 ±0.1) x 
10 45 erg s -1 , we expect Bq ~ 750 G at 10r g from the divergence theorem (Fig. [T5l) . 

For 0.97 < k r < 1.19, we get 1.1 < r core (43 GHz) < 0.01 pc and 0.18 < B x < 0.21 pc (us- 
ing Q ru at 43 GHz). An interesting feature of this source is that the estimated equipartition 
radius for an accretion-disk-launched jet (R eq ~ 0.5 pc) is within the range of VLBI observa- 
tions at 43 GHz assuming k r ~ 1. Multi-epoch VLBI monitoring of this jet at 43 GHz should 
provide evidence for the extended acceleration predicted from the magnetically driven jet 
model of K07 since the jet flow is not expected to reach its terminal Lorentz factor until a 
distance of ~ 1 pc for a black-hole-launched jet and ~ 5 pc for an accretion-disk-launched 
jet. 



3.4 0954+658 

In Fig. [16] we show the spectral index maps for 0954±658 from 4.6-43 GHz. We observe 
a strongly inverted spectrum from 4.6-22 GHz in the core region (Fig. [T6h . b, c), which 
then flattens considerably at the highest frequencies (Fig. [TBTd). We note the spectral index 



gradient across the jet in the 4.6-7.9 GHz image, which is discussed in 



O'Sullivan fc Gabuzda 
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Table 9. Averaged core-shifts using 43 GHz as reference frequency for 1418+546, plotted in Fig, 1131 



V 

(GHz) 


Ar 

(mas) 


projected 

(pc) 


Fraction of beam 
(%) 


4.6 


0.29 ±0.02 


0.75 ±0.05 


10 


5.1 


0.27 ±0.02 


0.70 ±0.05 


10 


7.9 


0.16 ±0.02 


0.42 ± 0.05 


9 


8.9 


0.13 ±0.02 


0.34 ±0.05 


9 


12.9 


0.085 ± 0.007 


0.22 ±0.02 


8 


15.4 


0.066 ± 0.007 


0.17 ±0.02 


8 


22.2 


0.036 ± 0.005 


0.09 ±0.01 


6 



Col. (1): observed frequency in GHz, Col. (2): magnitude of the shift from the 43-GHz core position, Col. (3): projected shift 
in parsecs (1 mas = 2.6 pc), Col. (4): core-shift as a fraction of the beam-size at the corresponding frequency, in percent. 

Table 10. Derived physical quantities for 1418±546 (fc r = 1) 





-Score {y) 


fcore {y") 


Bipc 




(GHz) 


(G) 


(PC) 


(G) 


(pc ■ GHz) 


4.6 


0.040 ±0.011 


4.1 ±2.3 


0.16 ±0.14 




5.1 


0.044 ±0.012 


3.7±2.1 


0.16 ±0.14 


3.7±2.0 


7.9 


0.067 ± 0.004 


2.7 ±0.3 


0.18 ±0.03 


4.1 ±0.4 


8.9 


0.076 ± 0.004 


2.2 ±0.2 


0.17 ±0.02 


3.8 ±0.4 


12.9 


0.11 ±0.01 


1.6 ±0.3 


0.17 ±0.06 


3.9 ±0.8 


15.4 


0.13 ±0.01 


1.5 ±0.2 


0.19 ±0.03 


4.4 ±0.5 


22.2 


0.19 ±0.02 


0.9 ±0.2 


0.17 ±0.04 


4.0 ±0.7 


43.1 


0.36 ± 0.02 


0.52 ±0.05 


0.19 ±0.03 


4.3 ±0.4 



Col. (1): observed frequency in GHz, Col. (2): radio core magnetic field strength in Gauss, Col. (3): de-projected distance of 
the radio core to the base of the jet in parsecs, Col. (4): magnetic field strength at 1 pc, Col. (5): average core-position offset 
measure, in units of 10 9 pcGHz, for shifts with vi in Table[8]equal to v in Col. (1) above. Note: Physical parameters at 4.6 GHz 
calculated with 5.1 GHz core-position offset measure. 



( 120091 ) . The region of inverted spectrum at ~3 mas from the core in the 8.9-12.9 GHz map 
(Fig. fT6b) may also be connected with the change in the jet direction in the plane of the 
sky. The extended emission in this source falls off sharply with increasing frequency, which 
is why we could not obtain core-shift measurements at frequencies higher than 8.9 GHz. 

We obtained two reliable core-shift measurements for the jet of 0954+658, with the angle 
of the shifts (~ 127°) consistent with the inner jet direction at 43 GHz of ~ —59° (Table 11). 
In this case, we assume a value of k r = 1 from 4.6-43 GHz due to our lack of measurements, 
and calculate the shifts with which to align our images from this. The spectral index maps 
in Fig. dn] indicate that this is a reasonable assumption between 4.6 and 43 GHz, providing 
a flat to inverted core spectrum with a steep spectrum jet as expected. 

In orde r to estimate the mag netic field strength f or this source, we use a Doppler factor 



of 6.2 from 



Hovatta et al. 



(I2009h and /3 app = 6.2 from 



Gabuzda et al. 



(I1994h to find T = 6.3 



and 9 = 9.3°. We get = 2.2° from the size of jet components within 1 mas of the core. 
We use a jet spectal index of -0.6 throughout the calculations due to the steep fall-off in the 
jet spectral index. In this case, to calculate the magnetic field strength we need to assume 



a value for the lower cut-off in the electron energy distribution since B\ oc 7 min • Using 
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Table 11. Measured core-shifts for 0954+658 







Ar 


Angle 


fraction of beam 


(GHz) 


(GHz) 


(mas) 


(°) 


(96) 


4.6 


7.9 


0.16 ±0.02 


126.8 ±0.6 


6 


4.6 


8.9 


0.18 ±0.02 


127.3 ±0.6 


7 



Col. (1) & (2): observed frequency in GHz, Col. (3): magnitude of the core-shift in milliarcseconds, Col. (4): angle of core-shift 
in degrees. Col. (5): core-shift as a fraction of the 4.6-GHz beam-size, in percent. 



Table 12. Derived physical quantities for 0954+658 



V 


Bcore (f) 


fcore {y) 






(GHz) 


(G) 


(pc) 


(G) 


(pc ■ GHz) 


7.9 


0.038 ± 0.002 


7.0 ±0.6 


0.27 ±0.03 


8.9 ± 0.7 


8.9 


0.044 ± 0.002 


6.0 ±0.4 


0.27 ±0.03 


8.7 ±0.6 



Col. (1): observed frequency in GHz, Col. (2): radio core magnetic field strength in Gauss, Col. (3): de-projected distance of 
the radio core to the base of the jet in parsecs, Col. (4): magnetic field strength at 1 pc, Col. (5): average core-position offset 
measure, in units of 10 9 pc ■ GHz, for shifts with V2 in Table [TT] equal to v in Col. (1) above. Note: magnetic field strengths 
calculated for 7 m i n = 100 



Vt TV = 8.8 ± 0.7 pc ■ GHz, we get Bi = 0.27 ± 0.03 G with 7 min = 100 while if 7 min extends 
all the way down to 1, we get Bi = 0.34 ± 0.05 G. Our estimates of the physical parameters 
for this source are shown in Table [121 assuming 7 m j n = 100 and a = —0.6. 

Extrapolating these measurements, we estimate a magnetic field strength in the 43 GHz 
radio core of B mrP . ~ 0.21 G at a distance of r core ~ 1.3 pc. Using a black hole mass of 10 8 ' 37 



M flWu et al J 12009( 1. we find R eq , A D ~ 0.2 pc and R eg (BH) ~ 0.02 pc. From the calculated 
total luminosity of (1.1 ± 0.2) x 10 46 erg s _1 , we expect B ~ 4 x 10 3 from the divergence 
theorem. Using B(r) = 0.27r _1 , we find B ~ 2 x 10 3 G at 10r r 



3.5 1156+295 and 1749+096 

For these two sources, we were unable to obtain core-shift measurements due to the weak 
nature of the jet emission. Hence, in this case we cannot be sure if the core corresponds to a 
stationary feature or the r = 1 surface. To produce the spectral index maps (Figs. [T7l[T8|) . we 
aligned the images by their model-fitted core positions. This produced uniform flat spectra 
in the core regions (e.g., the slice along the 4.6-7.9 GHz spectral index map for 1156+295, 
Fig.fTTh). We can see that the core spectrum of 1156+295, is strongly inverted (o^Jghz ~ 1-3) 
at low frequencies and then flattens considerably at the higher frequencies with o^ghz ~ 
—0.05 (Fig. [T7|) . The core spectrum of 1749+096 (Fig. [TBI displays similar properties, with 
^tJghz ~ 1.0 and becomes flat at the highest frequencies with «1|ghz ~ 0.01. 



4 DISCUSSION 



Magnetic Field Strength in AGN Jets 17 



In principle, the position of the vertex of the cone should provide a good estimate for the 
location of the central supermassive black hole powering the jet given that the distance to 
the radio core is likely to be much greater than the distance to the jet injection point or jet 
"nozzle" . However, the exact position is uncertain due to the strong dependence of r core on 
small deviations in the value of k r . In the case of 2200+420, the base of the jet can have 
a projected separation of up to r pro j. max ~ 0.4 mas (~ 0.5 pc for k r — 0.92) from the 43- 
GHz core. For 2007+777 and 1418+546, we find r projimax ~ 0.2 mas (k r = 0.91) and ~ 0.1 
mas (k r = 0.97), respectively, from their 43-GHz cores. Only the errors in the core-shift 
measurements are taken into account for our magnetic field strength calculations, hence, the 
results may be subject to systematic errors if one or more of the jet parameters ((f), 6, T, 5) 
are inaccurate. Of course all values presented here should be considered in the context of 
the assumption that the jet speed measurements accurately reflect the bulk flow speed of 
the jet and not some pattern speed. 

The main result from our core-shift measurements is that we find equipartition be- 
tween the particle and magnetic field energy densities to be an accurate assumption for 
the radio core of the observed sources from 4.6-43 GHz. Conservation of particle number 
in a freely expanding conical jet requires iV oc r~ 2 , which then implies B oc r _1 from 



magnetic field strength to 



equipartition. Furthermore, we expect the longitudina' 
as r~ 2 and the transver se component to fall-off as r _1 (IBegelman. Blandford fc Rees 



ail-off 



1984 



Hutter fc Mufson 



19861 ) so that over most of the length of the jet the transverse component 



will dominate and the overall magnetic field stren gth will scale as r l . Th is is also seen in 



models of magnetically dominated relativistic jets ( IVlahakis fc Koniglll2004l ) where, after the 
poloidal component of the magnetic field initially dominates, the toroidal component domi- 
nates at large distances and follows an approximate r _1 fall-off in field strength. Our results 
are well described by B cx r _1 , and by extrapolating back to the jet launching region, we get 
magnetic field strengths consistent with those required by models of magnetically launched 
jets to produce the observed jet luminosities. 

The implication of finding N oc r~ 2 consistent with our results is that, in the radio 
core from 4.6-43 GHz, continuous re-acceleration of the emitting particles is required to off- 
set both the radiation and adiabatic expansion losses in this region. This could be achieved 
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Table 13. Particle number density estimates at 1 pc for each source (k r = 1). 



Source 


ATi(7min = 1) 






(cm -3 ) 


(cm -3 ) 


0954+658 


1061 


6.8 


1156+295 






1418+546 


223 


2.2 


1749+096 






2007+777 


354 


3.5 


2200+420 


103 


1.0 



Col. (1): IAU source name, Col. (2): particle number density at 1 pc in cm 3 for a pair plasma with 7 m i n = 1, Col. (3): particle 
number density at 1 pc in cm -3 for a normal plasma with 7 m i n = 100. 



through numerous shocks (IMarscher fc Gearlll985l ). magnetic reconnection (IZenitani &: Hoshino 



20081 ) and/or turbulence from Kelvin- Helmholtz or current-driven instabilities. 

In order to decouple the magnetic field strength from the particle number density, we 
used the equipartition condition of equation [3j We can now put the calculated B\ back into 
the equation to get an estimate for the particle number density at 1 pc. However, we also 
need to assume a value for the lower cut-off of the electron energy distribution (7 m i n ). For 
a jet dominated by a no rmal (electron-proton) plasma, the energy distribution must have 



a cut-off of 7 min ^ 100 (jGhisellini et al. 



1992 



Ward 



c et al. 



19981), while 



1996 



or an electron- 



Hirotani 



20051 ) 



positiron jet 7 m i n can extend all the way down to 1 (iReynolds et al. 
In Table [TBI we list the particle number density in the jet at 1 parsec (iVi) for 7 m ; n equal to 
1 and to 100. 

Extended jet acceleration regions, most likely caused by MHD and/or hydrodynamical 



forces, have 



jeen detected in several AGN jets with the high resolution of VLBI (e.g., 



Jorstad et al. 



20051 ). From the model of K07, we can calculate the distance at which the jet 
approaches its terminal Lorentz factor. From this, we find that 2200+420, 2007+777 and 
1418+546 would be interesting candidates for the detection of extended jet acceleration from 
VLBI monitoring, since the terminal Lorentz factor distance for accretion-disk-launched jets 
for thes e sources extends be yond the distance of the 43-GHz VLBI core from the base of 
the jet. IJorstad et al.l (120051 ) have already detected accelerating components from 43 GHz 



VLBA monitoring of 2200+420 providing more weight to the estimates presented in this 
paper. 

Finally, we stress the importance of continued VLBI jet monitoring to obtain estimates of 
jet speeds, viewing angles and Doppler factors to enable magnetic field strength calculations 
for a much larger sample of AGN. 
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5 CONCLUSIONS 

From our measurements of the frequency dependent position of the VLBI radio core, we find 
that: 

1. The jets of 2200+420, 2007+777 and 1418+546 from 4.6-43 GHz are well described 
by the conical jet model of BK79 with r corc (z/) oc v~ x l kr . Furthermore, power-law fits to the 
data strongly support the equipartition value of k r — 1 in this region. 

2. Our spectral index maps are also consistent with the BK79 jet model of an optically 
thick, unresolved jet base with all maps displaying a flat or inverted spectrum in the core 
region, turning over to an optically-thin steep-spectrum region further downstream in the 
jet. 

3. We calculate equipartition magnetic field strengths with typical values from 10's to 
100's of mG in the 4.6-43 GHz radio cores. Extrapolation of our results, following an r _1 
power law, to the accretion disk and black hole jet-launching distances gives magnetic field 
strengths in general agreement with those expected from theoretical models of magnetically 
powered jets. 
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Figure 1. Spectral index maps for 2200+420 constructed from frequency pairs, with plots of slices along the jet region indicated 
by the solid line, (a) 4.6-7.9 GHz. Spectral index scale range: —1.2 to 1.5. Clip levels of / maps: 3.5 mjy (4.6 GHz), 3.8 mjy 
(7.9 GHz). Beam FWHM and Position Angle: 2.17x1.60 mas, -14° (4.6 GHz). Contours: 4.6 GHz /, starting at 3.48 mjy and 
increasing by factors of 2. (b) 8.9—12.9 GHz. Spectral index scale range: —1.2 to 0.7. Clip levels of / maps: 2.2 mjy (8.9 GHz), 
6.8 mjy (12.9 GHz). Beam FWHM and Position Angle: 1.21x0.84 mas, -14° (8.9 GHz). Contours: 8.9 GHz /, starting at 2.21 
mjy and increasing by factors of 2. (c) 15.4—22 GHz. Spectral index scale range: —2.3 to 1.0. Clip levels of / maps: 5.2 mjy 
(15.4 GHz), 5.9 mjy (22 GHz). Beam FWHM and Position Angle: 0.71x0.49 mas, -14° (15.4 GHz). Contours: 15.4 GHz /, 
starting at 4.35 mjy and increasing by factors of 2. (d) 22-43 GHz. Spectral index scale range: —0.8 to 1.0. Clip levels of / 
maps: 6.9 mjy (22 GHz), 13.3 mjy (43 GHz). Beam FWHM and Position Angle: 0.55x0.40 mas, -19° (22 GHz). Contours: 
22 GHz /, starting at 6.08 mjy and increasing by factors of 2. The arrows connect the plots of the spectral index slices to their 
corresponding maps. The hollow circle on the slice plot indicates the core position at the lower frequency of the pair. 
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Figure 3. Plot of core-shift measurements versus frequency for 2200+420 (listed in Table |3}, using 43 GHz as the reference 
frequency. Best-fit line: Ar = A^ -1 /*^ - 43. l" 1 /^) w i tn parameters A = 2.5 ± 0.3 and k r = 0.99 ± 0.07 
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Figure 5. Extended plot of B versus r for 2200+420 to include the extrapolated magnetic field strength {B ~ 0.14r~ 1 ) at the 
accretion disk (10r s ) and black holes (r g ) distances. The estimated equipartition radii from K07 are indicated by the arrows 
with -R cqi AD ~ 0.2 pc and i?. cq ,BH ~ 0.03 pc. Also plotted are points at 10r 3 (grey) and r g (black) representing the theoretically 
expected magnetic field strengths for a magnetically powered jet. 
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Figure 6. Spectral index maps for 2007+777 constructed from frequency pairs, with plots of slices along the jet region indicated 
by the solid line, (a) 4.6-7.9 GHz. Spectral index scale range: —1.1 to 1.0. Clip levels of / maps: 1.1 mjy (4.6 GHz), 2.0 mjy 
(7.9 GHz). Beam FWHM and Position Angle: 1.78x1.60 mas, 86.5° (4.6 GHz). Contours: 4.6 GHz /, starting at 1.07 mjy and 
increasing by factors of 2. (b) 8.9—12.9 GHz. Spectral index scale range: —1.4 to 1.1. Clip levels of / maps: 1.7 mjy (8.9 GHz), 
2.3 mjy (12.9 GHz). Beam FWHM and Position Angle: 0.95x0.84 mas, 82.5° (8.9 GHz). Contours: 8.9 GHz /, starting at 1.70 
mjy and increasing by factors of 2. (c) 15.4-43 GHz. Spectral index scale range: —1 to 0.6. Clip levels of / maps: 1.8 mjy (15.4 
GHz), 7.7 mjy (43 GHz). Beam FWHM and Position Angle: 0.55x0.50 mas, 83.6° (15.4 GHz). Contours: 15.4 GHz /, starting 
at 1.82 mjy and increasing by factors of 2. The arrows connect the plots of the spectral index slices to their corresponding 
maps. The hollow circle on the slice plot indicates the core position at the lower frequency of the pair. 
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Figure 7. Straight-line fit to the all the core-shift data for 2007+777, assuming equipartition. The y-intercept value is equal 
to —0.001 ± 0.004 mas. A straight-line fit through the origin is expected for k r = 1. 



0.25 




v [GHz] 

Figure 8. Plot of core-shift measurements versus frequency for 2007+777 (listed in Table [6}, using 43 GHz as the reference 
frequency. Best-fit line: Ar = A(i> _1 /*r - 43. 1" 1 /*"-) with parameters A = 1.08 ± 0.11 and k r = 0.97 ± 0.06 



Magnetic Field Strength in AGN Jets 27 




10 4 



100 



o 



0.01 



\I^(AD) 

V 



10" J 10" 4 °- 001 °- 01 ai 1 10 

r[pc] 

Figure 10. Extended plot of B versus r for 2007+777 to include the extrapolated magnetic field strength (B ~ 0.26r _1 ) at 
the accretion disk (10r s ) and black holes (r g ) distances. The estimated equipartition radii from the model K07 are indicated 
by the arrows with -R e q,AD ~ 0.7 pc and -R eqj BH ~ 0-1 P c - Also plotted are points at 10r 3 (grey) and r g (black) representing 
the theoretically expected magnetic field strengths for a magnetically powered jet. 
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Figure 11. Spectral index maps for 1418+546 constructed from frequency pairs, with plots of slices along the jet region 
indicated by the solid line, (a) 4.6—7.9 GHz. Spectral index scale range: —0.5 to 0.8. Clip levels of / maps: 1.7 mjy (4.6 GHz), 
2.6 mjy (7.9 GHz). Beam FWHM and Position Angle: 2.10x1.84 mas, 73.6° (4.6 GHz). Contours: 4.6 GHz /, starting at 1.66 
mjy and increasing by factors of 2. (b) 8.9—12.9 GHz. Spectral index scale range: —2 to 1.2. Clip levels of / maps: 2.1 mjy (8.9 
GHz), 3.7 mjy (12.9 GHz). Beam FWHM and Position Angle: 1.06x0.94 mas, 79.6° (8.9 GHz). Contours: 8.9 GHz /, starting 
at 2.14 mjy and increasing by factors of 2. (c) 15.4—22 GHz. Spectral index scale range: —1.0 to 1.0. Clip levels of I maps: 3.0 
mjy (15.4 GHz), 5.8 mjy (22 GHz). Beam FWHM and Position Angle: 0.69x0.54 mas, 61.9° (15.4 GHz). Contours: 15.4 GHz 
/, starting at 1.88 mjy and increasing by factors of 2. (d) 22-43 GHz. Spectral index scale range: —2.0 to 0.5. Clip levels of / 
maps: 4.2 mjy (22 GHz), 9.2 mjy (43 GHz). Beam FWHM and Position Angle: 0.48x0.42 mas, 80.8° (22 GHz). Contours: 22 
GHz /, starting at 4.14 mjy and increasing by factors of 2. The arrows connect the plots of the spectral index slices to their 
corresponding maps. The hollow circle on the slice plot indicates the core position at the lower frequency of the pair. 
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Figure 13. Plot of core-shift measurements versus frequency for 1418+546(listed in Table 1X31 , using 43 GHz as the reference 
frequency. Best-fit line: Ar = A(y~ x l h ^ - 43. l" 1 /^) witn parameters A = 1.4 ± 0.2 and k r = 1.08 ± 0.11 
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Figure 14. Plot of magnetic field strength (B) in units of Gauss versus distance along the jet (r) in units of parsecs for the 
core of 1418+546 at 4.6, 5.1, 7.9, 8.9, 12.9, 15.4, 22.2 Sz 43 GHz using k r = 1. Also plotted is the best-fit line B = Bir" 1 with 
B 1 = 0.178 G with a formal error of 0.006 G. 
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Figure 15. Extended plot of B versus r for 1418+546 to include the extrapolated magnetic field strength (B ~ 0.18r _1 ) at 
the accretion disk (10r s ) and black holes (r g ) distances. The estimated equipartition radii from the model K07 are indicated 
by the arrows with -R e q,AD ~ 0.5 pc and -R eqj BH ~ 0-1 P c - Also plotted are points at 10r 3 (grey) and r g (black) representing 
the theoretically expected magnetic field strengths for a magnetically powered jet. 
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Figure 16. Spectral index maps for 0954+658 constructed from frequency pairs, with plots of slices along the jet region 
indicated by the solid line, (a) 4.6—7.9 GHz. Spectral index scale range: —1.5 to 1.5. Clip levels of / maps: 1.5 mjy (4.6 GHz), 
3.8 mjy (7.9 GHz). Beam FWHM and Position Angle: 1.87x1.76 mas, -82.4° (4.6 GHz). Contours: 4.6 GHz I, starting at 
1.53 mjy and increasing by factors of 2. (b) 8.9—12.9 GHz. Spectral index scale range: —1.5 to 2.0. Clip levels of / maps: 2.3 
mjy (8.9 GHz), 2.6 mjy (12.9 GHz). Beam FWHM and Position Angle: 0.98x0.95 mas, 49.8° (8.9 GHz). Contours: 8.9 GHz 
7, starting at 2.26 mjy and increasing by factors of 2. (c) 15.4-22 GHz. Spectral index scale range: —2.0 to 2.0. Clip levels 
of I maps: 2.8 mjy (15.4 GHz), 7.8 mjy (22 GHz). Beam FWHM and Position Angle: 0.57x0.54 mas, -45.5° (15.4 GHz). 
Contours: 15.4 GHz /, starting at 2.31 mjy and increasing by factors of 2. (d) 22-43 GHz. Spectral index scale range: —1.0 to 
0.6. Clip levels of I maps: 8.9 mjy (22 GHz), 16.0 mjy (43 GHz). Beam FWHM and Position Angle: 0.37x0.36 mas, -12.9° 
(22 GHz). Contours: 22 GHz I, starting at 8.94 mjy and increasing by factors of 2. The arrows connect the plots of the spectral 
index slices to their corresponding maps. The hollow circle on the slice plot indicates the core position at the lower frequency 
of the nair. 
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Figure 17. Spectral index maps for 1156+295 constructed from frequency pairs, with a plot of a slice along the 4.6-7.9 GHz 
map indicated by the solid line, (a) 4.6-7.9 GHz. Spectral index scale range: —1.5 to 2.0. Clip levels of / maps: 2.7 mjy (4.6 
GHz), 10.9 mjy (7.9 GHz). Beam FWHM and Position Angle: 2.65x1.56 mas, 3.5° (4.6 GHz). Contours: 4.6 GHz I, starting 
at 2.73 mjy and increasing by factors of 2. (b) 8.9—12.9 GHz. Spectral index scale range: —1.0 to 1.0. Clip levels of / maps: 7.3 
mjy (8.9 GHz), 14.9 mjy (12.9 GHz). Beam FWHM and Position Angle: 1.43x0.83 mas, 1.1° (8.9 GHz). Contours: 8.9 GHz 
/, starting at 7.29 mjy and increasing by factors of 2. (c) 15.4-22 GHz. Spectral index scale range: —1.0 to 1.0. Clip levels 
of / maps: 5.8 mjy (15.4 GHz), 24.9 mjy (22 GHz). Beam FWHM and Position Angle: 0.80x0.47 mas, -0.3° (15.4 GHz). 
Contours: 15.4 GHz /, starting at 6.37 mjy and increasing by factors of 2. (d) 22-43 GHz. Spectral index scale range: —1.0 to 
1.0. Clip levels of / maps: 15.5 mjy (22 GHz), 44.2 mjy (43 GHz). Beam FWHM and Position Angle: 0.56x0.34 mas, 4.4° (22 
GHz). Contours: 22 GHz /, starting at 15.45 mjy and increasing by factors of 2. 
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Figure 18. Spectral index maps for 1749+096 constructed from frequency pairs, (a) 4.6—7.9 GHz. Spectral index scale range: 
-1.5 to 1.5. Clip levels of I maps: 2.4 mjy (4.6 GHz), 8.4 mjy (7.9 GHz). Beam FWHM and Position Angle: 3.18x1.46 mas, 
—0.3° (4.6 GHz). Contours: 4.6 GHz /, starting at 2.35 mjy and increasing by factors of 2. (b) 8.9-12.9 GHz. Spectral index 
scale range: -1.0 to 0.5. Clip levels of I maps: 6.4 mjy (8.9 GHz), 9.1 mjy (12.9 GHz). Beam FWHM and Position Angle: 
1.69x0.78 mas, -2.4° (8.9 GHz). Contours: 8.9 GHz I, starting at 6.41 mjy and increasing by factors of 2. (c) 15.4-22 GHz. 
Spectral index scale range: —1.0 to 0.5. Clip levels of I maps: 8.3 mjy (15.4 GHz), 15.7 mjy (22 GHz). Beam FWHM and 
Position Angle: 1.00x0.46 mas, —2.2° (15.4 GHz). Contours: 15.4 GHz J, starting at 5.95 mjy and increasing by factors of 2. 
(d) 22-43 GHz. Spectral index scale range: -1.0 to 0.5. Clip levels of I maps: 11.5 mjy (22 GHz), 39.3 mjy (43 GHz). Beam 
FWHM and Position Angle: 0.70x0.31 mas, -3.3° (22 GHz). Contours: 22 GHz I, starting at 11.54 mjy and increasing by 
factors of 2. 



